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Abstract: Changes in ecological stoichiometry reflect nitrogen (N), phosphorus (P) and both N and P 
limitations in a plant community, which in turn affect plant diversity of the community. However, the 
relationship between plant community diversity and ecological stoichiometry has not yet been fully 
researched in arid and semi-arid regions. Ecological stoichiometry and plant community diversity indices 
of eighteen communities in the upper reaches of Tarim River, northwestern China, were analyzed by 
multivariate analysis of variance in 2016. The correlation between ecological stoichiometry and plant 
community diversity was assessed by redundancy analysis (RDA). Results indicated that the Margalef 
index was significantly correlated with carbon (C) and P concentrations, the Simpson index and 
Shannon—Weaner index were significantly correlated with plant C concentration, and the Pielou index was 
significantly correlated with plant C and N concentrations. Moreover, C:N and C:P ratios had significant 
impacts on plant community diversity. Our results highlight the importance of ecological stoichiometry in 
driving plant community diversity in the upper reaches of Tarim River, northwestern China. 
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1 Introduction 


Ecological stoichiometry is an important tool for exploring the dynamic balance of various 
nutrient elements (Sterner and Elser, 2002; Zechmeister-Boltenstern et al., 2015; Moorthi et al., 
2016; Shang et al., 2018). At present, ecological stoichiometry has been applied to the molecules, 
cells, individuals, communities and ecosystems (Elser et al., 2007; Elser et al., 2010; Yan et al., 
2016). It has also been used to estimate the nutrient limitation of plants and to reveal the nutrient 
utilization strategies of plants (Sterner and Elser, 2002). The lack of carbon (C), nitrogen (N) and 
phosphorus (P) nutrients can limit plant growth (Sperfeld et al., 2017; Zeng et al., 2017), thus 
decreasing plant community diversity (Conti and Diaz, 2013). Plant community diversity is 
important for maintaining the stability and process of an ecosystem (Tilman et al., 1996; Naeem 
and Li, 1997; Yachi and Loreau, 1999; Snyder et al., 2006; Haddad et al., 2011; Pelini et al., 2014; 
Vilela et al., 2016; Zhang et al., 2017; Bustos-Segura et al., 2017; Saitta et al., 2018). Therefore, 
using ecological stoichiometry method to study plant communities and their relationships can 
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help to reveal the ecological strategies of plants in specific environments, and is of great 
significance to the material circulation and energy balance in a community. 

Beginning in the 1980s, ecological stoichiometry has become an important method for plant 
ecology research (Frost et al., 2002; Allen and Gillooly, 2009; Yan et al., 2016; Yu et al., 2017a; 
Castellanos et al., 2018; Li et al., 2019). Researches mostly focus on the effect of ecological 
stoichiometry on plant community diversity (Kerkhoff et al., 2005; Liu et al., 2018; Huang et al., 
2018) and examine the ecological stoichiometry of plant communities with different functional 
groups and changes in plant community diversity in grassland, forest and wetland ecosystems 
(Sardans et al., 2011; Yan et al., 2015; Tao et al., 2016; Soufbaf et al., 2017). However, little 
research was carried out on the desert ecosystem (Conti and Díaz, 2013). In addition, the study of 
ecological stoichiometry was always carried out at individual level not at the community level. 
(Halvorson et al., 2017; Zeng et al., 2017) 

In northwestern China, because climate change and anthropogenic disturbances occur 
frequently (Chen, 2013; Zhang et al., 2018), plant communities of this region are extremely 
vulnerable and are expected to suffer sharp reductions in plant community density in the future 
(Liu et al., 2005). Plant communities in the upper reaches of Tarim River located in an especially 
fragile ecological zone in northwestern China (Gong et al., 2016), yet not much is known about 
these plant communities. Hence, exploring the impact of ecological stoichiometry on plant 
community diversity provides a scientific basis for guiding vegetation restoration of the desert 
ecosystem in the upper reaches of Tarim River. Here, we hypothesis that (1) characteristics of 
ecological stoichiometry will affect the distribution of plants, thereby affecting plant community 
diversity of the upper reaches of Tarim River; and (2) distance from the Tarim River will affect 
the characteristics of ecological stoichiometry and plant community diversity indices, but 
relationship between ecological stoichiometric and plant community diversity was not affected. 


2 Materials and methods 


2.1 Study area 


The study area is located in the upper reaches of the Tarim River, northwestern China 
(40?25'-41?40'N, 80?10'-84?36'E; Fig. 1). The annual mean temperature is 10.4?C and the mean 
annual precipitation and evapotranspiration are 50.4 and 1800.0 mm, respectively (Gong et al., 
2016). The dominant plant species in this desert ecosystem are Populus euphratica, Halostachys 
caspica, Tamarix chinensis, Phragmites communis, Lycium ruthenicum and Glycyrrhiza uralensis. 


2.2 Methods 


Data were collected at three sampling zones in the upper reaches of Tarim River in August 2016 
(Fig. 1). The three zones were selected at the southern side of the river and named T1, T2 and T3, 
respectively. Within each sampling zone, six plant communities (i.e., sub-plots) were established 
to sample vegetation and the sub-plots were spaced at 1 km intervals so as to avoid problems of 
pseudoreplication. Thus, the total plant communities were eighteen and were coded (Table 1). 

Different sub-plot sizes were used to sample trees, shrubs and herbs. Trees were sampled in a 
20 mx20 m plot, while shrubs and herbs were sampled in 10 mx10 m and 1 mxl m plots, 
respectively. Approximately 50 g of leaf material from these individual plant species were 
collected and placed in envelopes with desiccant. After removing litter and stones, the top soil 
(0-15 cm) of these plant species was taken to laboratory and soil characteristics were analyzed. 
The vegetation composition and soil characteristics of different communities are shown in Table 
1. 


2.3 Leaf analyses 


Leaf samples were carefully cleaned and oven-dried at 60°C for 2 h. Soil samples were air-dried 
after being sieved using a 2-mm mesh size, and visible roots and organic debris were removed by 
hand. All samples were then ground to a fine powder for element analysis (Soil Agrochemical 
Analysis (Bao, 2000)). Total C content was determined by using the potassium dichromate 
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Fig. 1 Location of sampling zones in the upper reaches of Tarim River, northwestern China. T1, T2 and T3 are 
the sampling zones with six plant communities. 


Table 1 "Vegetation composition and soil characteristics of different communities 


Community Donsnantssécies Total Fre Cov Soil C Soil N SoilP Soil water 
code P number (g/kg) (g/kg) (g/kg) content (96) 


Populus euphratica, 


T11 : 168 1.39 0.69 6.58 0.59 0.63 12.57 
Halostachys caspica 

T12 P. euphratica, Tamarix chinensis 524 1.77 0.31 4.75 0.52 0.57 13.11 

yig A eiphratica, Lycium meen: zey L77 038 699 0.56 053 1423 


Glycyrrhiza uralensis 
T14 T. chinensis 363 1.33 0.31 6.87 0.54 0.61 14.19 


P. euphratica, T. chinensis, 


T15 : 109 1.33 0.38 1.19 0.22 0.58 12.86 
G. uralensis 
T16 P. euphratica 110 0.67 0.23 0.89 0.20 0.65 12.08 
T. chinensis, 
T21 Phragmites conatus 42 1.16 0.03 0.54 0.16 0.66 9.77 
T22 eM 36 L00 003 104 020 052 1176 
Hexinia polydichotoma 
Alhagi sparsifolia, 
T23 ad ; 81 0.55 0.28 0.56 0.22 0.45 9.51 
Karelinia caspia 
T24 T. chinensis, A sparsifolia 72 1.17 0.36 1.75 0.22 0.52 11.84 
T25 E euphratica; T. chinensis, 219 L72 028 204 019 052 1166 
A. sparsifolia, K. caspia 
T26 A.sparsifolia, 43 056 0.04 182 022 0.50 9.70 
P. euphratica 
Halostachys caspica, 
T31 Kalidium foliatum 80 1.11 0.63 4.66 0.38 0.50 6.24 
T32 H. caspica 8 1.05 0.25 7.39 0.59 0.56 8.38 
T33 T. chinensis, H. polydichotoma 45 1.00 0.25 7.33 0.68 0.49 9.19 
T34 Ho caspicd, A. sparsijoha, 336 161 028 4.28 0.62 0.51 5.05 
K. caspia 
T35 Tehimensis: A. sparsiyotia; 30 139 007 327 019 054 8.61 
K. caspia 
T36 T. chinensis, K. caspia 10 1.17 0.12 2.66 0.26 0.53 8.19 


Note: Fre, vegetation frequency; Cov, vegetation coverage; C, carbon; N, nitrogen, P, phosphorus. 
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external heating method (Bao, 2000), total N content was determined by using the Kjeldahl 
method (Bao, 2000) and total P content was determined by using the molybdenum-antimony 
colorimetric method (Bao, 2000). All leaf analyses were repeated three times. 


2.4 Alpha diversity indices 


Next, the alpha diversity indices were calculated by the following equations: 


D-1-YP', (D 
H'=-PinP, (2) 
Ma=(S-1)/InN, (3) 
E-H'InS, (4) 


where D is the Simpson index; H' is the Shannon-Wiener index; Ma is the Margalef index; E is 
the Pielou index; P; is the probability of occurrence of the it species in a community; S is the total 
number of species in a community; and N is the total number of individuals in a community. 


2.5 Ecological stoichiometry 


Leaf C, N and P concentrations were used to calculate the importance values of each community 
in trees, shrubs and herbs. According to these importance values, the stoichiometric eigenvalues 
for all of communities were weighted (Yan et al., 2008). 


2.6 Statistical analyses 


One-way analysis of variance (ANOVA) with a Duncan post hoc multiple comparisons test (Hsu, 
1996.) was used to compare the differences in C:N:P ratios among different communities. All 
statistical tests were assessed at P«0.05 level. Linear regressions were then performed to describe 
the relationships between C, N, P, C:N ratio, C:P ratio and N:P ratio. These analyses were 
conducted in SPSS version 19.0 software (Chicago, IL, USA). Redundancy analysis (RDA) was 
used to explore the correlation between plant community diversity index and C, N and P 
concentrations and their ratios. RDA was conducted in CANOCO 4.5 software (Wright, 2014). A 
Monte-Carlo test was performed to obtain the sequence of importance of variables. 


3 Results 


3.1 Characteristics of ecological stoichiometry in different communities 


The highest C concentration occurred in T13 community and the lowest in T31 community. The 
highest N concentration occurred in T23 community and the lowest in T14 community. The 
highest P concentration occurred in T34 community and the lowest in T22 community. All of 
highest stoichiometric ratios (i.e., C:N, C:P and N:P ratios) occurred in T1 zone and the lowest in 
T3 zone (Table 2). Significant differences in C, N and P concentrations and their ratios were 
observed within a zone and among different zones. In addition, N and P concentrations were not 
significantly different among T15, T25 and T35 communities (Table 2). 

Figure 2 showed the correlations among different C, N and P concentrations and their ratios. A 
strongly correlation was found between C and N concentrations, and between C and P 
concentrations (Fig. 2). A significant negative correlation was found between C:N ratio and N 
concentration, and between C:N ratio and N:P ratio. A strong and positive correlation was found 
between N:P ratio and N concentration and between C:P ratio and N:P ratio. However, an opposite 
pattern was found between C:P ratio and P concentration. 


3.2 Characteristics of plant community diversity index in different communities 


There were significant differences in the Simpson's index between the communities in T1 and T3 
zones. However, there was no significant difference in the other diversity indices between the 
communities in T1 and T3 zones (Fig. 3). For the Shannon-Wiener index, significant differences 
were observed between the communities in T1 and T2 zones, and in T1 and T3 zones, while no 
significant differences were observed between the communities in T2 and T3 zones (Fig. 3). For 
the Margalef index and the Pielou index, no significant differences were found among the 
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Table 2 Characteristics of C, N and P concentrations and their ratios in different communities 


M Er dum md C:N ratio C:P ratio N:P ratio 
T11 428.33+1.04?° 16.74+0.06“4 1.3240.07^* 26.05+1.948% 326.04419.39°° 12.53+0.3484 
T12 447.70X11.76^ — 18.58*0.96^*^ 1.1020.04^^ 24.1230.62P* 408.5149.94^^ 16.940.585 
T13 453.3943 .43^ 19.5940.92^* 1.09+0.024° 23.18x1.28** 414.3425.02^5^ 17.92x1.16^* 
T14 448.11+7.47^ 10.42+0.854 1.30+0.064* 43.1832.86^* 344.59+9.98 BP 7.9930.33** 
T15 434.0152.33^'  17.47+1.01^™ 0.67+0.09®° 24.95+1.30A%¢ 653.89+89.28"* 26.1442.3548" 
T16 449.0245.48^ . 16.8541.234°¢ 0.76::0.059* 26.7831.71P^ 596.76339.12^* 22.39+42.594° 
T21 451.3931.00^* 11.4021.17P^ 0.95:0.019* 39.8724.01^* 474.6234.679 11.99+1.308" 
T22 420.51+1.05° 14.643-1.04P^ 0.5940.01°%  28.82+1.20^P*° 714.55317.64^* 24.8932.26^ 
T23 441.61x10.20^* — 20.6742.26*" 0.92+0.10®° 21.5142.0254 486.68::59.06^* 22.93+4.94% 
T24 445.3549.884* 13.13+1.644° 1.07+0.06®° 34.24+4.078 418.7432.724 12.3441.65^^ 
T25 435.52x13.37^^ — 18.7632.57^* 0.69+0.02484 . 23,5623.81^" 629.96::34.05^B^ 27.0743.1148 
T26 400.73310.039 = 14.5141.92^P^ 1.36+0.03* 27.97+4.21 Beet 295.74413.95%* 10.70+1.338" 
T31 391.0541.46°° 11.35+0.84% 0.69+0.01°° 34.58+2.66" 562.7634.96^* 16.07+0.734" 
T32 405.3745.95?* — 11.62+2.00°" 0.71+0.02°° 35.6746.84^ 567.44+16.64"* 16.0520.72P* 
T33 432.]19*11.339 — 14.0741.595* 1.26+0.114° 35.22+4.884 345.24+24.51®° 11.31+2.23®° 
T34 433.15+22.08^° . 14.58x0.02^* 1.46+0.104* 29.72x].54P*» 297.79+6.26© 10.05+0.72®° 
T35 419.32411.914 — 17.93+2.8148 0.83+0.084° 23.80+43.934° 507.29x52.38P^ 21.5022.34 
T36 443.30418.88^' — 12.84+0.458" 1.3740.024” 34.5340.52^ 323.74419.57%° 9.3740.49** 


Note: Different uppercase letters within a column are significant differences among C, N and P concentrations and their ratios at 
different zones (T1, T2 and T3) in the same community at P«0.05 level. Different lowercase letters within a column indicate significant 


differences among C, N and P concentrations and their ratios at different communities in the same zone at P«0.05 level. 
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Fig. 2 Correlations among C, N and P concentrations and their ratios 
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Fig.3 Simpson's index (D; a), Shannon-Wiener index (A7; b), Margalef index (Ma; c) and Pielou index (E; d) in 
different plant communities. Different uppercase letters indicate significant differences among plant community 
diversity indices in T1, T2 and T3 zones at P«0.05 level. Bars indicate standard errors. T11-T16, T21—T26 and 
T31-T36 mean the six plant communities in T1, T2 and T3 zones, respectively. 


communities in T1, T2 and T3 zones. Overall, for all plant community diversity indices, higher 
values occurred in the plant communities in T1 zone than in T2 and T3 zones. 


3.3 Correlations between ecological stoichiometry and plant community diversity index 


RDA was used to examine the correlation between plant community diversity and ecological 
stoichiometry (Fig. 4). The first axis explained 64% of the variation, while the second axis only 


3 : 
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Fig. 4 Biplot diagram of redundancy analysis (RDA) for ecological stoichiometry and plant community 
diversity indices. In this biplot diagram, longer lines indicate stronger associations, and shorter lines indicate 
weaker associations. When the angle is between 0° and 90°, there is a positive correlation between two variables; 
when the angle is between 90° and 180°, there is a negative correlation between two variables; when the angle is 
90°, no correlation is observed between two variables. C, carbon; N, nitrogen; P, phosphorus; N, nitrogen; D, 
Simpson index; H', Shannon- Wiener index; Ma, Margalef index; E, Pielou index. 
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explained 2%. Only axis 1 reflected the relationship between ecological stoichiometry and plant 
community diversity index. The results of RDA indicated that C:N ratio and N concentration 
played important roles in explaining plant community diversity among all concentrations and 
ratios. The Simpson index was strongly and positively correlated with P concentration, and was 
weakly and positively correlation with C and N concentrations (Fig. 4). In contrast, C:N and N:P 
ratios had strong negative correlations with Simpson index, and C:P ratio showed a weak negative 
correlation (Fig. 4). A clear positive correlation was found between the Shannon-Wiener index 
and C, N and P concentrations (Fig. 4). Likewise, a negative correlation among these ratios was 
found, and the remaining indices showed similar patterns as well (Fig. 4). 

The effects of concentrations and ratios on plant community diversity index were different, and 
it is not possible to conclude which concentration or ratio had the greatest impact on plant 
community diversity index. Therefore, A Monte-Carlo test was performed on C, N and P 
concentrations and three ratios and obtained the sequence of importance of these stoichiometric 
variables. The importance sequence indicated that the most important element was N 
concentration (Table 3). 


Table 3 Sequence of importance and Duncan's test of carbon (C), nitrogen (N) and phosphorous (P) 
concentrations and their ratios 


Parameter Sequence Explanation F P 
N 1 0.624 14.632 0.002 
C:N ratio 2 0.571 13.302 0.002 
C:P ratio 3 0.539 12.583 0.002 
P 4 0.498 11.941 0.002 
C 5 0.258 4.532 0.041 
N:P ratio 6 0.103 1.524 0.235 


Figure 5 showed t-value results from RDA for the single stoichiometrical factor (i.e., N 
concentration, P concentration, C:N and C:P) influencing plant community diversity indices. The 
Simpson index, the Pielou index and the Shannon-Wiener index all fell within the circle and the 
Margalef index partly fell into the circle with solid line (Fig. 5a). This implied that N 
concentration strongly and positively correlated with the Pielou index, the Simpson index and the 
Shannon-Wiener index, and weakly and positively correlated with the Margalef index. In other 
words, with an increase in N concentration, the values of plant community diversity indices 
increased as well. The C:N ratio revealed similar patterns to those of N concentration (Fig. 5c). 
The Simpson index and the Shannon-Wiener index fell completely into the circle with solid line, 
and the Pielou index and the Margalef index partly fell into the circle with solid line (Fig. 5b). 
This result suggested that P concentration strongly and positively correlated with the Simpson 
index, the Pielou index and the Shannon-Wiener index, and weakly and positively correlated with 
the Pielou index and the Margalef index. And C:P ratio revealed similar patterns to those of C:N 
concentrations (Fig. 5d). 


4 Discussion 


4.1 Characteristics of ecological stoichiometry 


This study demonstrates large differences in C, N and P concentrations among different plant 
communities. C concentrations were highest at T13 community. This community has the largest 
number of Phragmites spp. that is a unique xerophilous plant (Tuboi and Hussain, 2018). Because 
T1 zone is the closest to the water source in the three sampled zones, a large number of 
Phragmites spp. grows in there. With the rapid growth of this species, the vascular bundles 
increase in the leaves and photosynthesis is enhanced, thereby leading to the effective 
accumulation of organic matter in the leaves (Behera et al., 2012). The highest N concentrations 
were observed in the species of Alhagi sparsifolia Shap. that dominates T23 community. A. 
sparsifolia is a naturally-growing and drought-tolerant plant that depends on groundwater (Xiang 
et al., 2018). If the plant is close to water source, its growth might be inhibited (Zhang et al., 
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Fig. 5 t-value results from RDA for single stoichiometrical factor influencing plant community diversity indices. 
(a), N concentration; (b), P concentration; (c), C:N ratio; (d), C:P ratio. The solid circle indicates positive 
correlation and the dotted circle indicates negative correlation. Plant community diversity indices fall into the 
solid circle, indicating a positive correlation with the elemental concentration or ratio. Plant community diversity 
indices fall into the dotted circle indicating a negative correlation with the elemental concentration or ratio. C, 
carbon; N, nitrogen; P, phosphorus; N, nitrogen; D, Simpson index; H', Shannon- Wiener index; Ma, Margalef 
index; £, Pielou index. 


2018). On the contrary, if the plant is far away from water source, it is likely that less water 
cannot meet its growth need. As such, the moderate amount of water is most beneficial for the 
growth of the plant. Moreover, A. sparsifolia is a leguminous plant that has high N fixing ability 
and utilization efficiency (Lei et al., 2014). All in all, the characteristics of ecological 
stoichiometry reflect the vegetation component of the community and the nutrient absorption and 
utilization of the vegetation. 

In general, N and P are assumed to be the two most limiting elements for plant growth (Elser et 
al., 2000; Elser et al., 2007; Qu et al., 2017). We found a significant difference in N:P ratios 
among different plant communities. The N:P ratio are higher than 16 in some communities, 
indicating that these communities were strongly P-limited (Yan et al., 2015; Yu et al., 2017b), and 
other communities that the N:P ratio are less than 14, were more limited by N (Plach et al., 2017). 
However, factors affecting plant N:P ratio are complex and comprehensive (Elser et al., 2007). 
The nutrient restriction in different communities are controlled by many factors, and the 
interaction between plant elements makes N:P ratio more complicated (Yan et al., 2015; Plach et 
al., 2017). Hence, the mechanism behind it needs further research. 


4.2 Characteristics of plant community diversity 


This study showed that the communities in T1 zone had the highest Simpson index values, 
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indicating that they had strong anti-interference abilities. Because of the intermittent erosion by 
the Tarim River for a long time, the communities are more resistant to interference (Sures et al., 
1999; Bu et al., 2017). However, other communities with poor anti-interference abilities have 
higher self-healing abilities after habitat destruction. The results indicate that dominant species 
can play a leading role in maintaining community structure and ecological function (Zhang et al., 
2017). The values of the Shannon-Wiener index of the T11 and T16 communities in T1 zone were 
high, indicating that the structure and composition of these communities are more complicated. In 
contrast, the T31 community in T3 zone had a low value of the Shannon-Wiener index, pointing 
to a simpler community structure. The T11 and T16 communities in T1 zone are located closer to 
the Tarim River and are not subjected to drought stress. The water resources are relatively 
abundant and thus more plant species could survive in there (Merritt et al., 2010). In contrast, the 
T31 community in T3 zone is far from the river, and only a handful of desert species that 
overcome drought stress could survive in there. The values of Margalef index and Pielou index in 
the plant communities of T1 zone were higher than those of the other zones, reflecting the high 
richness of species and abundant distribution of species in T1 zone. This result may be attributed 
to the fact that the plant communities in T1 zone are closest to the water source and have plenty of 
water that supports more species to survive and grow in there (Merritt et al., 2010). 


4.3 Correlations between ecological stoichiometry and plant community diversity index 


The results supported our hypothesis that ecological stoichiometry and plant community diversity 
were related. It is possible that the ecological stoichiometry is important in driving plant 
community diversity. Previous studies have found that C is the main element for plant, and 
compared with other nutrient, its concentration is high and stable (Heyburn et al., 2017), but C is 
not a key element that limits plant growth (Makino et al., 2003; Heyburn et al., 2007). The result 
showed no significant correlation between plant community diversity and C concentration was 
found, further confirming this latter view. In our study, N and P concentrations were important to 
plant community diversity indices that showed a significant correlation with N and P 
concentrations. Study has showed that changes in N and P concentrations had an impact on plant 
community diversity (Zhao et al., 2010). Both the increase and decrease of N and P 
concentrations will affect the growth rate of plants, thus changing plant community diversity 
(Frenken et al, 2017). The sequence of importance of ecological stoichiometry on plant 
community diversity shows that C:N and C:P ratios also had significant impacts except for N 
concentration. The C:N and C:P ratios provide insight into plant's ability to absorb nutrients and 
assimilate C (Bell et al., 2018), to a certain extent, they can also reflect the nutritional efficiency 
of plant (Liu et al., 2010). Plant communities have different C sequestration efficiencies, while C 
accumulation efficiency and storage capacities are associated with N and P supplies that limit 
plant growth. Ultimately, ecological stoichiometry and plant community diversity appeared to be 
closely linked (Castellanos et al., 2018). 


5 Conclusions 


The importance of ecological stoichiometry in driving plant community diversity within a desert 
ecosystem was confirmed in this study. Firstly, the closer to the riparian, the greater the 
stoichiometric concentration of the community was, especially the C and N concentrations. 
Secondly, the closer to the riparian, the greater the plant community diversity index was. However, 
the plant community diversity index showed an opposite pattern when the distance to the riparian 
continues increase. What's more, RDA showed that N concentration was a key factor affecting 
plant community diversity in the Tarim River. In conclusion, under extreme environment of 
lacking nutrients and water, desert plants in Tarim River have formed special characteristics of 
ecological stoichiometry to maintain plant community diversity. 
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